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■Our  study  considered  effective -plane -wave  propagation,  both 
longitudinal  and  shear,  through  a  medium  containing  a  random  distribution  of 
spherical  inclusions.  We'.  assumed  that  the  particles  and  matrix  are 
separated  by  a  thin  layer \ of  elastic  material  with  different  properties. 

For  some  systems,  we  predict  measurable  effects  for  the  thin  layers. 
Especially,  we  considered  P^/epoxy  and  (ilC/fidP. 


1.  Introduction 

Wave  propagation  through  a  partible -reinforced  composite  medium  has 
been  studied  by  many  authors  [1-13].  Except  [5],  all  these  studies  assume 
that  the  Inclusions  bond  perfectly  with  the  surrounding  matrix  material.  In 
[5],  for  long  wavelengths,  the  authors  consider  the  effect  of  a  thin  viscous 
third  layer.  Recently,  Sayers  [11]  examined  the  effect  of  this  layer  when 
the  particles  and  the  matrix  possess  the  same  properties. 

In  the  present  study,  we  analyze  the  problem  of  wave  propagation  in  a 
composite  medium  with  a  random  distribution  of  spherical  inclusions.  The 
inclusions  are  separated  from  the  matrix  by  thin  layers  of  elastic  material. 
The  properties  of  the  layers  vary  through  the  thickness  such  that  there  is  a 
continuous  transition  from  the  inclusions  to  the  matrix. 

The  object  of  this  study  was  to  explore  the  practicality  of  using 
ultrasound  to  characterize  properties  of  interface  layers.  Ultrasound  is  a 
practical  tool  for  measuring  properties  of,  and  characterizing  the  state  of, 
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a  material  with  microstructure  (or  changes  in  microstructure).  References 
to  such  studies  occur  in  ftirfX 


2.  Scattering  by  a  spherical  inclusion  with  an  interface  layer 

Consider  a  spherical  inclusion  of  radius  a  and  elastic  properties  Aj , 
Pi,  Pi  embedded  in  an  elastic  matrix  of  material  properties  A2 ,  p2>  Pi- 
Also,  let  the  inclusion  be  separated  from  the  matrix  by  a  thin  layer  of 
uniform  thickness  h(«a)  with  variable  material  properties  A(r),  p(t) ,  and 
p( r) .  Here,  A,  p  denote  Lam6  constants  and  p  density.  Let  A(r),  p(r)  be 
expressed  as 


A(r)  +  2p(r)  -  (Ax  +  2/i1)f(r),  a  <  r  <  a  +  h, 


(1) 


p(r)  -  Pig(r),  a  <  r  <  a  +  h.  (2) 

>• 

Here,  f(r)  and  g(r)  denote  general  functions  of  r.  A  special  case  arises 
when 


f (a)  -  1,  g(a)  -  1. 


£(* +  h)  -  •  *<* +  h>  *  S  • 


(3) 


with  the  stipulation  that  f(r)  and  g(r)  with  their  first  derivatives  are 
continuous  in  (a,  a  +  h) .  Since  h  is  assumed  to  be  much  smaller  than  a,  it 
follows  from  (3)  that  f'(a)  and  g’(a)  can  be  approximated  by 


f  •  (a)  _  C^2  +.tp2)  (A,  +  2/ii ) 

* }  h(Ax  +  2pi ) 
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Another  special  case  arises  when  the  interface  material  possesses 
constant  properties.  Then  we  have 


■X'c-S 

:\ 


fH 


□  □ 


3 


I 

v!' 

,( - 


!i* 

V 


iii1 

$ 

f* 

<>* 

L 

a 

A 


ivl 

|I7| 

**J| 


a 

1 


§ 

§ 


$ 

IS 


f (r)  -  <AJ  +  2mD/(A1  +  2mi),  g(r)  -  mI/Mi-  (5) 

Here  A£t  n[  are  the  Lam6  constants  for  the  interface  material. 

Ve  also  make  the  assumption  that  h  is  much  smaller  than  the  wavelength 
of  the  propagating  wave.  Then,  to  first  order  in  h/A,  A  being  the 
wavelength, 

t  s  ,  i  t  s  ,  i 

r  —  r  +  r  ,  r_  —  r.  +  r., 
rr  rr  rr  r©  r©  r© 


t  _  s  i 
rr$  ”  rr#  rr$ 


Here  rjj  is  the  stress  tensor  and  superscripts  t,  s,  and  i  denote  the 
transmitted,  scattered,  and  incident  field  quantities,  respectively.  Note 
that  rf j ,  r[j,  and  tJj  appearing  above  are  calculated  at  r  -  a.  The 
spherical  polar  coordinates  r,  8,  *  are  defined  in  Fig.  1.  Boundary 
conditions  (6)  express  the  fact  that,  to  first  order  in  h/A,  the  traction 
components  do  not  suffer  any  jump  across  the  layer.  However,  the 
displacement  components  suffer  jumps  given  by 
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Here , 


K  .  <*x  K  -  f1  _ — 
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Using  equations  (3)  and  (4)  in  (10),  we  find  that  Kj  and  K2  are  given 
approximately  by 
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On  the  other  hand,  if  eq.  (5)  is  used,  then 

Kx  -  (Ax  +  2mi)/(AJ  +  2M{).  K2  -  Mi/Ml-  (13) 

Mai  and  Bose  [5]  studied  a  problem  similar  to  the  one  considered  here.  They 
assumed  a  thin  viscous  third  layer  between  the  sphere  and  the  matrix  and 
imposed  the  condition  of  radial -displacement  continuity. 

We  assumed  the  incident  wave  to  be  either  a  plane  longitudinal  wave 
propagating  in  the  positive  z-direction  or  a  plane  shear  wave  polarized  in 
the  x-direction  and  propagating  in  the  positive  z-direction.  Thus, 
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Here,  ^  -  u/cx  and  k3  -  w/c3.  u  denotes  the  circular  frequency  of  the  wave 
and  Ci,  c3  denote  the  longitudinal  and  shear  wave  speeds  in  the  matrix.  The 
factor  e”lwt  was  suppressed. 

ul  given  above  can  be  expanded  in  spherical  vector  wave  functions  as 
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Vector  wave  functions  end  N^  appearing  above  are  regular  at  r  - 

0  and  are  given  by 
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The  scattered  and  transmitted  fields  can  be  written 
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Here,  the  prime  denotes  that  kt  and  k2  are  to  be  replaced  by  k}  (-  w/cj)  and 
k2  (-  w/c2) ,  respectively,  cj  and  c2  are  the  wave  speeds  in  the  inclusion.  . 

L^3\  M^3\  and  N^3^  are  obtained  by  replacing  by  in  (16).  Note  that 

Jn  is  the  spherical  Bessel  function  of  the  first  kind,  and  hR  is  the 

spherical  Hankel  function  of  the  first  kind. 

The  constants  A,  B,  C,  A’,  B' ,  C'  are  found  by  using  conditions  (7)  and 
(8) -(10).  For  this  purpose,  we  define  the  following  matrices: 


Fn(kxa)  -  nhn(kxa)  -  kxahn+1 (kxa) , 

Gn(k2a)  -  n(n  +  1)  hn(k2a) ,  Hn(kxa)  -  hn(kxa), 
Ixx(k2a)  -  (n  +  1)  hn(k2a)  -  k2ahn+i(k2a)  , 

SFn(kxa)  -  (n2  -  n  -  4  k?a2)  hn(kxa)  +  2kxahn+i (kxa) 
SGn(k2a)  -  n(n+l)  [ (n  -  1)  hn(k2a)  -  k2ahn+i (k2a) ] , 
SHn(kxa)  -  (n  -  1)  hn(kxa)  -  kxahn+i(kxa) , 

SIn(k2a)  -  (n2  -  1  -  hk?a2)  hn(k2a)  +  k2ahn+1 (k2a) . 

Equations  to  determine  Ann  and  are  found  to  be 
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M„,  are  obtained  from  Mn,  respectively,  by  replacing  hn  and  hn+j  by 
jn  and  Jn-t-i,  respectively,  in  (19)  and  (20),  and  by  replacing  kj  and  k2  by 
k{  and  k2 ,  respectively .  In  writing  (21)  we  express  u£  and  given  by  (15) 

as 
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It  then  follows  that 


i  , 

rr(mn) 

•  -  l  if-1 

r  1 

ur(mn) 

li 

a  n  n 

li 

'■rr0(mn)J 

>u0(mn)‘ 

(23) 


Ln  and  Mn  are  obtained  from  and  respectively,  by  replacing  hn  and  hn+i 
by  jn  and  jn+1 ,  respectively. 

The  equation  to  find  B^  is 
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Here , 
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Once  Anjn,  Cmn,  and  are  determined  by  solving  (21)  and  (24),  the 
scattered  field  is  then  found  from  (17).  Since  the  expressions  for  the 
field  inside  the  inclusion  will  not  be  needed  to  derive  the  dispersion 
equation  governing  the  effective  wavenumber  of  plane -wave  propagation 
through  the  composite  medium,  we  omit  these. 

3.  Distribution  of  inclusions 

In  [5,  14]  the  scattered-field  expressions  were  used  to  calculate 
effective  wave  speeds  at  long  wavelengths  in  a  medium  with  a  distribution  of 
spherical  inclusions  with  interface  layers.  A  'quasicrystalline' 
approximation  together  with  the  assumption  of  no  correlation  was  used  to 
derive  expressions  for  the  effective  wave  speeds.  As  has  been  shown  [12] , 
particular  forms  of  two-particle  correlations  can  be  included  in  the 
formalism.  But  this  leads  to  complicated  equations  that  require  numerical 
solutions . 

In  this  study,  we  adopt  the  approach  taken  in  [10,  11]  to  calculate 
approximate  phase  velocities  and  attenuation  of  plane  longitudinal  and 
plane-shear  waves.  In  this  simple  approximation,  the  effective  wavenumber  k 
is  related  to  the  forward  scattered  amplitude  by  the  equation  [15]: 

K2  -  kg  +  4irnf (K)  .  (26) 

Here,  k0  is  the  wavenumber  in  the  absence  of  scatterers,  n  the  number 

density  of  scatterers,  and  f  the  averaged  forward- scattered  amplitude. 
Equation  (26)  is  an  implicit  equation  for  the  determination  of  the  (complex) 
wavenumber  K.  A  further  simplification  occurs  when  the  solution  to  (26)  is 
taken  as 

K2  -  kg  +  4»rn0f(k0)  .  (27) 

Equation  (27)  was  derived  by  Foldy  [16]  and  has  been  used  by  many  authors  to 
calculate  the  frequency  dependence  of  phase  velocity  and  attenuation  of 
plane  waves.  Equations  (26)  and  (27)  are  valid  for  low  volume 
concentrations  of  inclusions. 


Using  (17)  in  (26),  we  find  that  for  longitudinal  waves  the  effective 
wave  number  is 
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For  shear  waves  we  obtain 
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In  the  following  section  we  present  phase  velocity  and  attenuation 
calculated  from  the  above  two  equations.  Note  that  (he  real  part  of  K/k 
gives  the  velocity  ratio  c/C,  where  c  is  the  velocity  in  the  matrix  and  C  is 
the  effective  velocity  in  the  composite.  The  attenuation  of  power  is 
obtained  from  the  equation 


g  -  2  Im  £  . 
k  k 


(30) 


For  dilute  concentration  this  equation  reduces  to 


g  -  n°  2- 


(31) 


Here,  E  denotes  scattering  cross  section. 

4.  Numerical  results  and  discussion 

Computations  were  made  for  two  particular  composite  materials:  lead- 
epoxy  and  SiC-Al.  For  both  materials,  we  consider  two  interface 
thicknesses:  zero  and  0.1a,  through  which  the  properties  vary  linearly  from 

the  inclusion  to  the  matrix. 

Numerical  results  for  the  lead-epoxy  composite  based  on  a  simplified 
eq .  (27)  were  presented  in  [10],  Figures  2  and  3  show  the  attenuation  and 


phase  velocity  of  a  longitudinal  wave  when  the  interface  thickness  equals 
zero.  These  results  agree  with  those  given  in  [10].  Also  shown  in  these 
figures  are  the  results  for  h/a  -  0.1.  It  is  seen  that  this  thickness  has 
little  effect.  At  intermediate  frequencies,  the  interface  lowers  the  phase 
velocity  and  slightly  increases  the  attenuation.  Figures  4  and  5  show  the 
results  for  the  shear  wave;  the  effect  is  slightly  larger  at  moderate 
frequencies.  We  also  computed  the  phase  velocity  and  attenuation  using  eq. 
(26).  Figure  6  shows  that  eq.  (27)  overestimates  the  attenuation.  For 
phase  velocity,  however,  as  shown  in  Fig.  7,  eq.  (26)  overestimates  it  at 
low  frequencies,  underestimates  it  at  moderate  frequencies,  and  gives  nearly 
the  same  results  at  high  frequencies. 

Finally,  in  Figures  8-11  we  show  results  for  the  second  example: 

SiC— Al.  These  results  are  based  on  eq.  (27)  and  show  that  the  interface 
decreases  both  the  attenuation  and  the  phase  velocity. 


5.  Conclusions 

We  considered  the  effect  of  thin  interface  layers  between  the 
inclusions  and  the  matrix  in  modifying  the  dynamic  properties  of  composite 
materials.  Dynamic  effective  properties  were  calculated  by  using  Foldy's 
equations.  We  found  that  interface  effects  are  larger  in  some  composites. 

We  also  studied  the  predictions  based  on  iterative  solutions  of  modified 
Foldy's  equations  in  which  the  scattering  amplitudes  were  calculated 
assuming  that  the  matrix  had  the  properties  of  the  composite.  This 
iterative  solution  underestimates  attenuation  in  general;  at  low  frequencies 
it  overestimates  phase  velocity. 
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Table  1. 

Properties  of  constituents 


Density  (g/cm3) 

E(GPa) 

M(G Pa) 

Lead 

11.3 

23.57 

8.35 

Epoxy 

1.18 

4.31 

1.57 

SiC 

3.181 

440.6 

188.1 

A1 

2.706 

71.6 

26.7 

Fig.  1.  Geometry  of  a  spherical  inclusion  with  an  interface  layer. 


Fig.  2.  Attenuation  of  longitudinal  waves  in  a  lead-epoxy  composite  with 
and  without  interface  layer. 

Fig.  3.  Phase  velocity  of  longitudinal  waves  in  a  lead-epoxy  composite 
with  and  without  interface  layer. 

Fig.  4.  Attenuation  of  shear  waves  in  a  lead-epoxy  composite  with  and 
without  interface  layer. 

Fig.  5.  Phase  velocity  of  shear  waves  in  a  lead-epoxy  composite  with  and 
without  interface  layer. 

Fig.  6.  Comparison  of  attenuation  coefficients  for  longitudinal  waves  in  a 
lead-epoxy  composite  predicted  by  an  iterative  solution  of  eq. 

(26)  and  by  eq.  (27). 

Fig.  7.  Comparison  of  phase  velocities  of  longitudinal  waves  in  a  lead- 

epoxy  composite  predicted  by  an  iterative  solution  of  eq.  (26)  and 
by  eq.  (27). 

Fig.  8.  Attenuation  of  longitudinal  waves  in  an  SiC-Al  composite  with  and 
without  interface  layer. 

Fig.  9.  Phase  velocity  of  longitudinal  waves  in  an  SiC-Al  composite  with 
and  without  interface  layer. 

Fig.  10.  Attenuation  of  shear  waves  in  an  SiC/Al  composite  with  and  without 
interface  layer. 

Fig.  11.  Phase  velocity  of  shear  waves  in  an  SiC/Al  composite  with  and 
without  interface  layer. 


iterative 


Lead-Epoxy 
h/a  ■  0.0 
c  *0.15 


Sayers,  Smith 
(1983) 


